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Abstract 

The linear instability of Gortler vortices in compressible boundary 
layers is considered* Using asymptotic methods in the high wavenumber regime, 
it is shown that a growth rate estimate can be found by solving a sequence of 
linear equations. The growth rate obtained in this way takes non-parallel 
effects into account and can be found much more easily than by ordinary 
differential equation eigenvalue calculations associated with parallel flow 
theories. 
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1. INTRODUCTION 


Our concern is with the linear growth of Taylor-Gortler vortices in com- 
pressible boundary layers. We develop a simple method for generating curves 
of constant amplification rates in the high wave number regime. 

The growth of Taylor-Gortler vortices in incompressible boundary layers 
has received a lot of attention in recent years due to its relevance to 
Laminar Flow Control (see for example Harvey and Pride [1]). The original 
calculation by Gortler [2] showed that Taylor's [3] instability mechanism 
which occurs for curved flows also operates in boundary layer flows. However, 
the relative complexity of the basic state for a boundary layer flow makes it 
a much more difficult task to examine the instability of this state theoreti- 
cally. Thus, the essential difficulty with the linear instability problem is 
that the growth of the boundary layer cannot in general be ignored and the 
appropriate linear instability equations are therefore partial differential 
equations. 

The original calculation by Gortler ignored the effect of boundary layer 
growth completely and his numerical results were later corrected by Hammerlin 
[4] who found that instability occurs first at zero wavenumber. Later calcu- 
lations by Hammerlin [5] and Smith [6] attempted to remedy this deficiency by 
including higher order curvature terms or terms associated with the non- 
parallel nature of the basic state. Further work by Herbert [7] for example 
was aimed at understanding why the various linear theories did not give con- 
sistent results. 

Floryan and Saric [8] gave a multiple scale approach to the linear 

• • 

Gortler instability problem along the lines of, for example, the work of 
Gaster [9] or Saric and Nayfeh [10] for Tollmien-Schlichting waves. Thus, 
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Floryan and Saric derived the partial differential equations governing the 
growth of Gortler vortices* These equations had been given in a more general 
context some years earlier by Gregory, Stuart and Walker [11] who discussed 
the instability of three-dimensional boundary layers. The equations can also 
be inferred from the work of Smith [6]. The solution given by Floryan and 
Saric [8] followed the approach of previous investigations and implicitly made 
a parallel flow approximation. By the latter phrase we mean that some intrin- 
sic property of the nonparallel nature of the basic state was ignored in solv- 
ing the disturbance equations. In fact, the above authors replace streamwise 
partial derivatives of the vortex by local spatial growth rates thus reducing 
the system to a set of ordinary differential equations. It is not clear how 
such an approach can be justified but when the growth rate vanishes the solu- 
tion can be interpreted as a local Taylor series solution of the full partial 
differential equations. The relevance of the solution elsewhere is not 
immediately apparent. 

More recently, Hall [12-13] has shown how asymptotic and numerical 
methods can be used to take non-parallel effects into account in a self- 
consistent manner. In the first paper, it was shown that small wavelength 
Gortler vortices can be described asymptotically using a multiple scale 
method. Hall found that the vortices locate themselves so as to maximize 
their local spatial growth rate. This requires that the vortices are con- 
centrated in a viscous layer in the interior of the flow. 

In the linear regime, Hall [13] solved numerically the full partial dif- 
ferential instability equations at 0(1) wavenumbers. The linear equations 
were found to be parabolic in the streamwise direction so that an initial dis- 
turbance was imposed at some location and its development followed as the 
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equations were marched downstream. The growth of the disturbance was followed 
by calculating the local rate of change of a disturbance energy density. The 
neutral position was defined to be the location where this local growth rate 
vanished. Not suprisingly, it was found that this position was a function of 
the location and form of the initial disturbance. Thus it was concluded that 
there exists no unique neutral curve for the Gortler problem. However, at 
high wavenumbers, the different neutral curves merge into the asymptotic and 
parallel flow neutral curves. The same would be true for the different 
possible growth rate curves. 

Thus, in the only regime where analytical progress is possible, the 
growth rate can be written down in asymptotic form and no numerical eigenvalue 
calculations are required. It is this idea which we will now apply to com- 
pressible boundary layers to show how growth rates for these flows can be 
simply calculated. 

Previous calculations of the compressible Gortler problem have used the 
parallel flow assumption to reduce the instability problem to an eigenvalue 
problem associated with an eighth order differential system. (See for example 
Aihara [14], Kobayashi and Kohama [15] or El-Hady and Verma [16].) In 
particular, El-Hady and Verma formulated the linear stability problem along 
the lines of Floryan and Saric [8] and gave curves of constant growth rate for 
various flow conditions. We show how these curves can be generated much more 
simply in the only regime where they are meaningful. The method we use is 
based on the asymptotic theory of Hall [12] for the incompressible problem. 
The method can be easily used for any flow configuration and needs little 
computational power. The method is based on the assumption that the vortex 
wavelength is small compared to the boundary layer thickness. The range of 


-4- 


validity of the methods can only be checked by a numerical solution of the 
full partial differential system governing the growth of vortices in growing 
compressible boundary layers. However, in general the high wavenumber regime 
is ultimately always applicable to any constant wavelength disturbance vortex 
developing in a growing boundary layer and so is therefore always physically 
relevant. 

The procedure adopted in the rest of this paper is as follows: In Sec- 
tion 2, we formulate the partial differential system governing small Gortler 
vortex disturbances in compressible boundary layers. In Section 3, we solve 
these equations for large wavenumbers and determine the spatial growth rates 
of the disturbances. Finally, in Section 4 we present our results and draw 
some conclusions. 


2. FORMULATION OF THE INSTABILITY EQUATIONS 

Apart from some minor differences, our formulation is essentially the 
same as that of El-Hady and Verma and so the reader is referred to that paper 
for more details. We choose L to be a typical streamwise length scale and 
take v , U , p , T , u to be the scales for the kinematic viscosity, 

00 00 OO 00 oo 

velocity, density, temperature, and coefficient of viscosity respectively. If 

* 1 x* 

the curvature of the wall at the streamwise location x is kCj-) * we 

define the curvature parameter 6 by 

r _ L 

6 " A * 

and a Reynolds number R by 

UL 

OO 

R = 

\> 

00 

•* 

and consider the limit R " with the Gortler number 


( 2 . 1 ) 

( 2 . 2 ) 



(2.3) 


G = 



6 


held fixed. The free stream Mach number M is defined by 

00 ' 


M 



(2.4) 


where y and R are the ratio of specific heats and gas constant re- 
spectively. We define (x,y,z) to be dimensionless variables in the stream- 

wise, normal and spanwise directions scaled on L, R - ^^L, and R“^^L 
respectively. We shall assume that the vortices grow spatially in the x 
direction and therefore we consider them to be steady. 

The basic flow is written in the form 


(u,v,w) 


U o (u(x,y), 


R' 1/2 v(x,y), 0) 


+ • • • 


T = T oo T u ■ y(x,y). 


— 2 — 

P = P 00 p(x,y), p = p oo U oo p(x) 

where 

P* (u u x + v u y } = (y u ) - P x (2.5a) 

(pu) + (pv) = 0, (2.5b) 

x y 

p T x + v T y ] = rly (i*y y + ( y-!> ^ W u y 2 , (2.5d) 

T = T w , u = v = 0, y = 0, 


(2.5e) 
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together with condition on T and iT as y ». Here r is the Prandtl 
number whilst T w is the wall temperature. 

We now perturb (2.5) to a disturbance periodic in the z direction. The 
velocity components of the disturbance and temperature are scaled in an 
identical manner to the corresponding basic state quantities. The linearized 
instability equations take the form 


£ (u U) + ya 2 U + £ vU - (yU ) + £ u V 


t y y y x y 


[£(uu + vu)+(yu)]T-yu T =0, 

t 2 x y y y y y 


(2.6a) 


£ (v + kuG)U - cy U - (c+1) yU - y U 
Y x y x xy p x y 


+ 2 - (vV) + P - V + y a 2 V - (c+2) (yV ) + P 

t y^ x yyy 


[— _ (u v + v v + 1/2 kGu 2 ) + (c+1) y u + c y u 
Y 2 x y xy y x 


+ (c+2) (yv) +yu]T-yu T 
yy xy y x 


[cy u x + ( c+2 ) y v y ] T y - cy y iaW - (c+l)iay W y = 0, (2.6b) 


y iaU + (c+1) yia U + y ia V + (c+1) yiaV - iaP 
x x y y 


+ cy (u + v ) iaT - — pW - (c+2) ya 2 W - — pW 

x y T x Y y 


+ W y)y ■ °» 


(2.6c) 
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< HrK + <?)„ + !• A - T 

ij x iji y ip 


- (£ =P\ ■ 0 


(2.6d) 


- T U - 2(y-l) M 2 y u U + — T V - [^t (u T + vT) 

Y x “ y y Y y t 2 x y 


+ (y-1) M 2 £ u 2 +I(Il)]T+-oT + £ a 2 T 


y r y y T x r 


+ (— v-^-yT)T - | (u T ) =0. 

Y r p y' y r yy 


(2.6e) 


Here y = - - - and a is the spanwise wavenumber. We note that the above 

d T 

equations can be simplified if p is independent of x; in that case, we 

can set p = 1 in which case (2.6) reduce to the equations of El-Hady and 

*• 

Verma. The Gortler number G is defined by (2.3), c = X/y where X is the 

bulk viscosity. Equations (2.6) are to be solved subject to the perturbation 

quantities vanishing at y = 0, «° . 


3. THE HIGH WAVENUMBER SOLUTION FOR M ~ 0(1) 

00 

It is known from the work of Hall [12-13] that small wavelength Gortler 
vortices are located in the boundary layer so as to maximize their local 
spatial amplification rate. For the incompressible case and zero amplifica- 
tion rate, this position corresponds to where Rayleigh's criterion is most 

— 1/2 

violated. The depth of this layer is 0(a ' ) so we define n by 


n = (y - y(x)> a 


(3.1) 
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where y(x) is the as yet undetermined location of the layer. The stream- 
wise disturbance velocity in this layer expands as 

— 1 /2 2 x 

u = [U Q (n,x) + a UjOijX) + •••] exp {a / a(x) dx} 

where 

a(x) = 8 Q (x) + a + ••• 

Similar expansions hold for Va -2 , Wa~^ 2 , Pa - "^ 2 , and T whilst Z. 
expands as 

u = Uq (x) + UjP a + ••• . 


Again, similar expansions for T, v and y hold. The details of the expan- 
sion procedure are essentially identical to those of Hall [12] so we shall 
omit a lot of detail here. The Gortler number G expands as 


G = g Q a . 


(3.2) 


It is convenient to define the matrix A by 


A(x,y) = 


- Vo 

Hn + 


0 T 
1 


g 0 U 0 


0 

0 


u. 


Mo + 


Vo 


0 T 


" g 0 U 0 

=2“ 

2 T o 

*0 . Vo 
r T o 


(3.3) 
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Here p has been set to unity since for the asymptotic solution, pressure 
may be rescaled by local edge pressure. If the above expressions are sub- 
stituted into the disturbance equations and like powers are equated, we obtain 
the system of equations 


A(x,y) 


U, 


0 


T J 

0 J 


= 0, 


(3.4) 


at zeroth order whilst the next two order systems yield 


U 

A(x, y) {V 


(U, 


3A , -v 

" _n 97 (x ’y } 


'01 


and 


("*) 2 - 1 

( U 0nn| 

I 1 

(°o| 

1 V 1 


' \7 l 

ui I 

in 

1 • -1- 1 

It") 


(3.5a) 


(3.5b) 


0 


The coefficient matrices C, E, and F can be written down in terms of 
quantities involving basic flow quantities. The system (3.4) has a nontrivial 
solution if 

|A| = 0 (3.6) 


which for a given choice of g^ and y determines three possible spatial 
amplification rates 3^. The first order solution can then be written as 


U C 

V, 


= u Q (n ,x) 


- T 0 r Vo + 

— I— + w 0 ] 

U 1 T o 


-o> 


(3.7) 


0 


- r u o e o 

U 1 [ ~=~ 
T 
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= u o - • 


The system 


* TR + n 
A Uq = 0 


will then have a solution and will be needed at higher order. Here, 


u 0 = 


, u o e o T o 

•{ + Un> 


‘0 J - 


l o 


u o g o 


- r 5 ° 6 ° + " \ 
- u 0 {— + U 0 } 

0 

- - , u o B o ^0, 

2 Vir- + -> 

T o 


At this stage y(x) and U q (ti ,x) remain undetermined but at next 
order we find that (3.5) has a solution only if 


/ + n TR 9A ( — v n 

(u 0 ) 9? (x ’ y) i = ° 


(3.8) 


and this fixes the location y(x). Physically (3.8) can be interpreted as 
the condition that (5 q has a maximum at the layer y - y. The solution 
of (3.8) can then be written in the form 


n U U, 

O 


where (a, 3) satisfies 
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!^ +( ^. + !^ )6 = ( 0 ,! i ! 

T o T o V 

V , 9 o Vi , - ,Vi 

T T 1 T 1 T 

A 0 0 0 l 0 


2T 2 8 0 "0 Y 

T T T 

0 0 0 


2u 2 )/T q , 0) • a 


V - ^ 


(3.9) 


- ^ n) (7) - T <n) (y) 

u = , T = , etc. 

n n! ’ n n! ’ 


Finally, Uq is determined when the required solvability condition is 


applied to (3.6); this yields 


e u 0 n„ + U 0 + x 6 l"o - °- 


(3.10) 


Here e, y, and X are defined by 


e = (uJ) TR E a , y = (uJ) TR C a, X = (u£) TR F a . (3.11a,b,c) 


The coefficient matrices E, C and F are given in Appendix A. The solu- 
tions of (3.10) which decay to zero when |q | « can be written down in 
terras of parabolic cylinder functions, the most unstable one is 


U Q - exp {-1/2 n 2 /-J4 


(3.12) 


and the corresponding eigenrelation is 


= /-Ye /X. 
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If we were interested in finding the neutral Gortler number, we would have 
included higher order terms in (3.2) and set = 0. In that case, 

3 

(3.13) would be replaced by an equation to find g^, the order a term 
in the expansion of the neutral Gortler number. 

We now summarize the steps required to find 8^ and 8^ — the first 

two terms in the expansion of the spatial amplification rate. Firstly, at any 
depth y, the cubic equation specified by the condition 

I A | =0 

is solved for the three possible values of 8 q« The value of y is then 

varied until (3.8) is satisfied and then 8^ is determined by (3.13). 

Thus, it is not necessary to solve any differential equations numerically to 
obtain Bq and 8^* The answer we obtain is formally valid when a » 1, 
at smaller values of a it is at least as valid as the solution of El-Hady 
and Verma which would require large amounts of computer time. At progressive- 
ly higher values of a, the different approaches will converge. At finite 
values of the wavenumber a full numerical solution of (2.6) along the lines of 
Hall [13] is required. We now turn to the results we have obtained using the 
above approach. 


4. RESULTS AND DISCUSSION 

In Figure 1, we have shown curves of equal spatial amplification rate at 
a Mach number of 2 in the wavenumber - Gortler number plane. These curves 
correspond to the adiabatic wall condition being applied to the temperature. 
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The stagnation temperature is held at 311°K. The results shown are formally 
valid at large values of a and are then the only unique amplification rates 
which exist for the Gortler pro blem . In this figure the wavenumber is made 

/v 

non-dimensional by using L = / . Similar curves for Mach number of 4 

e 

and 6 are given in Figures 2 and 3 respectively. We note that the constant 
growth rate curves extend to small wavenumbers for higher Mach numbers. How w 
ever, the small wavenumber region is perhaps beyond the range of validity of 
the theory. The neutral curves for all three Mach numbers are shown in Figure 
4. These curves shift towards the left with increasing Mach number, indicat- 
ing the stabilizing effect of compressibility. 

Next, we perform a calculation which is of direct engineering signifi- 
cance. We consider the boundary layer on the wall of a supersonic nozzle (see 
reference [17]). Transition in this boundary layer is caused by Gortler 
vortices. The flow accelerates to Mach 3.5 towards the exit of the nozzle* 
The distribution of the local edge Mach number, Hartree's pressure gradient 

parameter (g, ), and the Gortler number based upon momentum thickness 0 
h 

x 


are plotted in Figure 5. The amplification factor (N = f a(x)dx) for 
Gortler vortices is computed for a fixed physical wavelength (local nondiraen- 
sional wavenumber is also plotted in the figure) using the current asymptotic 
theory and the parallel theory used in the computations of reference [17]. 
For design purposes, the agreement between the two approaches is fairly good 
but it should be noted that the parallel theory calculations requires at least 
30 times as much computer time when compared with the asymptotic calculations. 
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APPENDIX A 

The non-zero elements of the coefficient matrices C, 
3.11) are 


T o T o 


8 0 ( - Vl -..h. 2 T 2- 1+ a , T l“l 

c n ■ -“2 - =- '”2 - -=- + u o a =-> - =- ,] + - 
L o o L o 0 

1 - Vl _ ? 1 2 Y 2 

c 19 = it f“ 3u T + 2 -zr 1 - Vt r> - zr>l 

T o T o T o T o 




u,T, 


T, 2 T„ 


5 0 r - . 11 - ,, K 

c 2i = — [_u 2 + ~ u 0 (( ~ } " ~ )] 

T o T o T o T o 

- B , + e JL ff - ^2), 

T o T o 

. 8 Vo Vi, 

+ —j— (u - — ) 

T T 

0 0 


0, 


u,T, 


T, 2 T. 


'0 r- “n . - ,,*lv" 2, , 

c 22 - -"2 - 7- [u 2 - + ' 7 -11 

L o 0 0 0 


-g 


'23 


'31 


2T' 


0 ro - 
2 
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u.T 


T. 2 T„ 


u t i l *1 ~ *2 _ ^1 u 0 

[2u n {u - J-l + u [(r±) - -i]} + (u - - Tr -: 

U <4 *nn 7p T 


T? 


A i T i u n 6 n u i T i 

. . i_ [2 T, - -Li - 6(<u, - -L2) ^ - -Ui 

T o T o T o T o 


3T T 

r _ ro 7 12 % 

C 32 = [3 3 ? 2 

*0 0 T 


Vl.-„ T U 2 T 2„ V 


«33--lT £+ ^~ + W -^ )} # ] 


T 0 T 0 


J 11 


E, F (Equation 


2 u 2 J 
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Figure 5. Variation of local Mach number (M e ), Hartree pressure gradient 
parameter ((3^), Gortler number based on momentum thickness (G^), wave 
number (a) and N factors in a Mach 3.5 supersonic nozzle. 
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